The escape of heavy atoms from the ionosphere of 
HD209458b. II. Interpretation of the observations 



T. T. Koskinen^'*, R. V. Yelle^, M. J. Harris'', P. Lavvas'^ 

Lunar and Planetary Laboratory, University of Arizona, 1629 E. University Blvd., 

Tucson, AZ 85721, USA 
^Department of Physics and Astronomy, University College London, Gower Street, 

London WCIE 6BT, UK 
"Groupe de Spectrometrie Moleculaire et Atmospherique UMR CNRS 6089, Universite 
Reims Champagne- Ardenne, 51687, France 



Abstract 

Transits in the H I 1216 A (Lyman a), O I 1334 A, C II 1335 A, and 
Si III 1206.5 A lines constrain the properties of the upper atmosphere of 
HD209458b. In addition to probing the temperature and density profiles 
in the thermosphere, they have implications for the properties of the lower 
atmosphere. Fits to the observations with a simple empirical model and a 
direct comparison with a more complex hydrodynamic model constrain the 
mean temperature and ionization state of the atmosphere, and imply that 
the optical depth of the extended thermosphere of the planet in the atomic 
resonance lines is significant. In particular, it is sufficient to explain the ob- 
served transit depths in the H I 1216 A line. The detection of O at high 
altitudes implies that the minimum mass loss rate from the planet is ap- 
proximately 6 X 10^ kg s~^. The mass loss rate based on our hydrodynamic 
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model is higher than this and imphes that diffusive separation is prevented 
for neutral species with a mass lower than about 130 amu by the escape of 
H. Heavy ions are transported to the upper atmosphere by Coulomb colli- 
sions with H+ and their presence does not provide as strong constraints on 
the mass loss rate as the detection of heavy neutral atoms. Models of the 
upper atmosphere with solar composition and heating based on the average 
solar X-ray and EUV flux agree broadly with the observations but tend to 
underestimate the transit depths in the O I, C II, and Si III lines. This 
suggests that the temperature and/or elemental abundances in the thermo- 
sphere may be higher than expected from such models. Observations of the 
escaping atmosphere can potentially be used to constrain the strength of the 
planetary magnetic field. We find that a magnetic moment of m < 0.04 mj, 
where mj is the Jovian magnetic moment, allows the ions to escape globally 
rather than only along open field lines. The detection of Si^"^ in the ther- 
mosphere indicates that clouds of forsterite and enstatite do not form in the 
lower atmosphere. This has implications for the temperature and dynamics 
of the atmosphere that also affect the interpretation of transit and secondary 

eclipse observations in the visible and infrared wavelengths. 
Keywords: Extra-solar planets, Aeronomy, Atmospheres, composition. 

Photochemistry 
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1 1. Introduction 



The detection of H, O, C"*", and Si^"*" in the upper atmosphere of HD209458b 



3 (Vidal-Madjar et al. 2003, 2004 Linsky et al. 2010), and the tentative de 



ifection of H in the upper atmosphere of HD189733b (Lecavelier des Etangs 



et al.l 2010) and Mg+ in the upper atmosphere of WASP-12b (Fossati et al. 



2010) are among the most exciting recent discoveries related to the atmo- 



spheres of extrasolar giant planets (EGPs). The observations demonstrate 
that the the upper atmospheres of close-in EGPs such as HD209458b differ 
significantly from the thermospheres of the giant planets in the solar system. 
They are much hotter, they extend to several planetary radii and instead of 
molecular hydrogen, they are primarily composed of atoms and atomic ions. 

The detection of heavy atoms and ions such as O, C+, Si2+, and Mg+ 
implies that the atmospheres of close-in EGPs are not always affected by 
diffusive separation. A likely explanation is that diffusive separation of the 
heavy atoms and ions is prevented by momentum transfer collisions with 
the rapidly escaping light atoms and ions. Mass fractionation during hy- 
drodynamic escape is believed to have played an important role in the early 



Involution of the atmospheres of the terrestrial planets (e.g., Zahnle and Kast 



ing, 1986; Hunten et al. , 1987) but it cannot be observed in action anywhere 



in the solar system. Observations of EGP atmospheres thus provide a unique 
opportunity to study this phenomenon that should lead to a better under- 
standing of evolutionary processes in our own solar system. 

Extended thermospheres give rise to large transit depths in UV trans- 



24 mission spectra. However, they are also potentially detectable in optical and 



25 infrared spectra. For instance, Coustenis et al. (1997, 1998) searched for an 



exosphere around 51 Peg b in the near-IR. In line with the current under- 
standing, they suggested that the exospheres of close-in EGPs such as 51 
Peg b are hot and composed primarily of atoms and ions. They also argued 
that hydrodynamic escape can lead to the escape of heavier species, and 
give rise to large in-transit absorption by such species in optical and near-IR 
spectra. However, 51 Peg b is not a transiting planet and the search for an 
exosphere around it was not successful. On the other hand, once the transit 



33 of HD209458b was first detected ( Charbonneau et al. , 2000 Henry et al 



2000), similar searches on this planet were also undertaken. 



Moutou et al. (2001) looked for visible absorption by species such as Na, 



36 H, He, CH+, C0+, N+, and H2O+ in the upper atmosphere of HD209458b. 



These observations were followed by Moutou et al. (2003) who attempted 



38 to measure the transit depth in the He 1083 nm line that was predicted 



39 to be significant by Seager and Sasselov (2000). The most recent searches 



40 were reported by Winn et al. (2004) and Narita et al. (2005) who looked for 



transits in the visible Na D, Li, Ha, H/5, H7, Fe, and Ca absorption lines. So 
far none of the ground-based searches have led to a detection of the upper 
atmosphere. However, the non-detection is based on only a few observations 
that have proven difficult to analyze, and the search should continue. 

Visible and infrared observations have been able to probe the atmosphere 
of HD209458b at lower altitudes. In fact, HD209458b was the first EGP 
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to have its atmosphere detected by transmission spectroscopy. The first 



48 detection was achieved by Charbonneau et aL (2002) who observed a deeper 



in-transit absorption in the Na D 589.3 nm resonance doublet compared to 
the adjacent wavelength bands. This detection was based on four transits 
observed with the Space Telescope Imaging Spectrograph (STIS) onboard the 



Hubble Space Telescope (HST) (Brown et al. , 2001). The same data were 



later reanalyzed by Sing et al. ( 2008a|[b ) who combined them with other 



observations (Knutson et al. , 2007) and created a transmission spectrum of 



HD209458b at wavelengths of 300-800 nm. They argued that the abundance 
of sodium in the atmosphere is depleted above the 3 mbar level either by 
condensation into Na2S clouds or ionization. We note that the detection of 



Si^^ in the thermosphere (Linsky et al. , 2010) constrains cloud formation 



59 mechanisms in the upper atmosphere and implies that the depletion of Na 



at low pressures is probably due to ionization (see Section 3.3) 



The atmosphere of HD209458b has also been observed several times in 



62 the infrared with the Spitzer space telescope. Deming et al. (2005) detected 



the secondary eclipse of the planet at 24 /im by using the Multiband Imaging 
Photometer (MIPS). Together with similar observations of TrES-1 at 4.5 and 



65 8.0 /im obtained by Charbonneau et al. (2005) who used the Infrared Array 



Camera (IRAC), these observations constitute the first detections of infrared 



67 emission from extrasolar planets. They were followed by Richardson et al 



68 (2007) who observed the infrared emission spectrum of HD209458b between 



69 7.5 and 13.2 fj,m with the Infrared Spectrograph (IRS). This spectrum was 
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reanalyzed by Swain et al. (2008) who noted that it is largely featureless with 



some evidence for an unidentified spectral feature between 7.5 and 8.5 /im. 



Knutson et al. (2008) used IRAC to observe the secondary eclipse of 



HD209458b in the 3.6, 4.5, 5.8, and 8.0 fim bands. They observed a higher 
than expected flux in the 4.5 and 5.8 fim bands and interpreted this as ev- 
idence for the presence of a stratospheric temperature inversion that gives 



76 rise to strong water emission at these wavelengths. Beaulieu et al. (2010) 



observed the transit of the planet in the same wavelength bands and also 
found evidence for the presence of water vapor in the atmosphere. The de- 
tection of water vapor is also supported by Swain et al. ( 2009[ ) who used the 
Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on HST 
to observe the secondary eclipse of HD209458b between 1.5 and 2.5 fim. In 
addition to water vapor, the NICMOS observations revealed the presence 
of methane and carbon dioxide in the emission spectrum. These detections 
provide valuable clues to the overall composition of the atmosphere but the 
uncertainties in the data and degeneracies between temperature and abun- 
dances in the forward model prevent a more quantitative characterization of 
the density and temperature profiles. 

In general, difficulties associated with reducing the data and the need 
to describe a few uncertain data points with models of growing complexity 
has led to disagreements on the analysis and interpretation of transmission 
and secondary eclipse data on different exoplanets. The same is true of the 



92 FUV transit observations of HD209458b. Vidal-Madjar et al. (2003) used the 
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STIS G140M medium resolution grating to detect a 15 ± 4 % transit depth 
in the wings of the stellar H Lyman a emission lin^ Based on the data, they 
argued that the planet is followed by a cometary tail of escaping hydrogen 
atoms that are accelerated to velocities in excess of 100 km s~^ by stellar 



ijradiation pressure (see also Schneider et al. 1998). Later, Vidal-Madjar et 



98 al. (2004) used the STIS G140L low resolution grating to detect absorption 



99 by H, O and C"*" in the upper atmosphere, and argued that the atmosphere of 



HD209458b escapes hydrodynamically. Linsky et al. (2010) used the Cosmic 



Origins Spectrograph (COS) on HST to confirm the detection of C^ and 
reported on the detection of Si^"*" around the planet. They also argued that 
the atmosphere escapes hydrodynamically. 



Ben-Jaffel (2007, 2008) disagreed with the interpretation of the G140M 



observations. He reanalyzed the G140M data and argued that the 15 % 
H Lyman a transit depth was exaggerated because the data were partly 
corrupted by short-term variability of the host star and geocoronal Lyman a 
emissions. He also showed that there is no evidence for a cometary tail in the 
data, and that absorption by H in the extended thermosphere of the planet 



no can explain the observations. Ben-Jaffel and Hosseini (2010) reanalyzed the 



[pl40L data and reached a similar conclusion regarding H (see also Koskinen 



et al. , 2010a). However, Ben-Jaffel and Hosseini (2010) fitted the H Lyman a 



113 observations by scaling the density profiles from the model of Garcia Munoz 



^The core of the H Lyman a line is entirely absorbed by the interstellar medium (ISM) 
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114 (2007), and argued that suprathermal O and C"*" are required to fit the transit 

115 depths in the O I and C II lines. They did not explain how the suprathermal 

116 atoms form and simply chose their properties to fit the observations. 



Holstrom et al. (2008) offered yet another explanation for the H Lyman a 



118 observations. They argued that hydrogen atoms cannot be sufficiently accel- 

119 erated by stellar radiation pressure before they are ionized by stellar X-rays 

120 and EUV (XUV) radiation. Instead, they suggested that the observed ab- 

121 sorption arises from a cloud of energetic neutral atoms (ENAs) that form by 

122 charge exchange between the protons in the stellar wind and the escaping 

123 hydrogen. In their model, the observed absorption reflects the velocity of 

124 the stellar wind, and the data can potentially be used to characterize the 



125 magnetosphere of the planet. Ekenback et al. (2010) recently updated the 

126 model to include a more realistic description of the escaping atmosphere and 

127 stellar wind properties. However, both studies ignored absorption by H in 

128 the thermosphere, which is significant, and did not attempt to explain the 

129 presence of heavier atoms and ions such as O, C^, and Si^^ in the escaping 

130 atmosphere. 

131 All of the interpretations of the UV transit data require that HD209458b 



132 has a hot and extended thermosphere. In fact, Koskinen et al. (2010a) (here- 

133 after KIO) also showed that absorption by thermal H and O in such a ther- 

134 mosphere explains both the H Lyman a and O I 1304 A transit depths. 

135 Further, their model related the observations to a few physically motivated 

136 characteristics such as the mean temperature and composition of the upper 
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137 atmosphere. They also used the observations to constrain the pressure level 

138 where H2 dissociates and estimated that the H2/H transition should occur 

139 at 0.1-1 /ibar. However, these results are based on an empirical model that 

140 was simply designed to fit the observations. One of the aims of the current 

141 paper is to show that the results are also supported by more complex physical 

142 models. 

143 We have also attempted to establish a more comprehensive description 

144 of the thermosphere that treats it as an integral part of the whole atmo- 

145 sphere rather than a separate entity. In order to do so, we developed a new 

146 one-dimensional escape model for the upper atmosphere of HD209458b that 

147 includes the photochemistry of heavy atoms and ions, and a more realistic de- 



ws ipcription of heating efficiencies. This model is described in detail by Koskinen 



et al. ( 2012[ ) (hereafter. Paper I). We also used results from a state-of-the-art 



150 photochemical model (Lavvas et al., in preparation) to constrain the density 

151 profiles of the observed species in the lower atmosphere. We discuss the im- 

152 plications of our results in the context of different observations, and show 

153 that observations of the upper atmosphere also yield valuable constraints on 

154 the properties of the lower atmosphere. 

155 2. Methods 

156 2.1. Stellar emission lines 

157 The interpretation of the FUV transit measurements relies on accurate 

158 characterization of the stellar emission lines and, provided that the ISM is 
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optically thick over parts of the line profile, their absorption by the ISM. The 
observed transit depths also depend on stellar variability. In this section we 
discuss the properties of the line profiles, stellar variability, and absorption 
by the ISM. The properties of the H Lyman a and the O I 1304 A triplet 
lines (hereafter, the O I lines) have been discussed in detail before (e.g.. 



Ben-Jaffel and Hosseini 2010 Koskinen et al. 2010a), and we present only 



a brief summary of them here. However, the COS observations of the C 
II 1335 A multiplet (hereafter, the C II lines), and the Si III 1206.5 A line 



(hereafter, the Si III line) (Linsky et al. , 2010) have not been modeled before. 



168 and thus we constructed emission line models for these lines based on the new 

169 data. These, and the previous models for the H Lyman a and O I lines, were 

170 used to calculate the predicted transit depths in Section [3j The observed 

171 line-integrated transit depths for IID209458b are listed in Table [2j 

172 2.1.1. H Lyman a and O 1 1304 ^ lines 



Wood et al. (2005) measured the H Lyman a emission line of HD209458 



174 by using the high resolution echelle E140M grating on STIS, and used the 

175 details of the line profile to constrain the column density of H in the ISM 

176 along the line of sight (LOS) to IID209458. Following KIO, we adopted the 



177 reconstructed line profile and ISM fit parameters from Wood et al. (2005) for 



178 the transit depth calculations in this paper. The short-term variability of the 



179 H Lyman a emissions from HD209458 was estimated by Ben-Jaffel (2007) 



ISO who used the G140M data to derive a magnitude of ~8.6 ± 5.6 % for this 
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variability - although the large uncertainty of the observations makes it diffi- 
cult to separate variability from noise. There are no observations to constrain 
the long-term variability or center-to-limb variations of the Lyman a line on 
HD209458. However, typical solar characteristics provide some guidance on 
these properties. 



Woods et al. (2000) studied the variability of solar Lyman a emissions 



based on satellite observations spanning four and a half solar cycles between 
1947 and 1999. They found that the variability ranges between 1 and 37 
% during one period of solar rotation (27 days), and the average variability 
during one solar rotation was found to be 9 ± 6 %. This result agrees well 
with the estimated variability of the Lyman a emissions from HD209458. The 



192 rotation period of HD209458 is estimated to be ~10-11 days (Silva-Valio 



2008). The G140M observations covered three different transits and took 



place within a month and a half. Each observation covered approximately 
2 hours in time. Thus the data can be affected by short-term variability 
and it is essential that such variability be properly accounted for. For this 



197 reason, we compare our models with the results of Ben-Jaffel (2007, 2008) 



198 and Ben-Jaffel and Hosseini (2010) who analyzed the data in the time tag 



mode and accounted for variability before calculating transit depths. 

Short-term variability in the H Lyman a line is mostly related to plage 
activity that is modulated by solar rotation while long-term variability de- 



pends on variations in both plage and active network regions (Woods et al. 



2000 ) . In fact, it is misleading to imagine the transits in any of the FUV emis- 
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204 sion lines as the planet crossing a smooth, uniformly emitting disk. Rather, 

205 one has to imagine the planet crossing a relatively dark disk with scattered 

206 bright regions. For instance, during solar maximum the plage^ cover ap- 

207 proximately 23 % of the solar disk while active network covers about 8.5 %. 

208 The brightness contrasts of the plages and active network are 6.7 and 3.2, 



209 Respectively, when compared to the quiet sun (Worden et al. , 1998 Woods 



210 et al. , 2000). This means that the transit depth in the H Lyman a can vary 

211 by factors of 0.4-3 as a function of time during maximum activity. Similar 

212 variability can be expected in the other FUV emission lines. 

213 The transit depth depends mostly on the path of the planet across the 

214 stellar disk although the plage and active network coverage can also change 

215 slightly with stellar rotation during the transit. This highlights the impor- 

216 tance of light curve analysis to work out transit depths - limited 'snapshots' 

217 during transit may be corrupted by the planet either covering or not cover- 

218 ing an active region. However, it is also possible for the transit depth to be 

219 altered even if the variations are not immediately evident in the light curves. 

220 This complicates the analysis of the observations, and underlines the need 

221 for repeated observations at different times. We note that the plage and 

222 active network coverage is significantly smaller during solar minimum than 

223 it is during solar maximum. This means that the transit depth is likely to 

224 be more stable and closer to the true transit depth during stellar minimum. 



^We include enhanced network in the definition of a plage here. 
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Unfortunately, the activity cycle of HD209458 has not been studied in detail, 
and thus we have no information about it. 

As we noted above, the properties of the O I lines were discussed exten- 



sively by KIO who fitted parameterized solar line profiles (Gladstone, 1992) 



229 to the O I lines of HD209458 that were observed with the STIS E140M grat- 



230 ing (Vidal-Madjar et al. , 2004). We adopted the O I line profile and ISM 



parameters from KIO for this study, and do not discuss these properties fur- 
ther here. We note that the plage and active network contrasts of the O 



233 I lines are 4.2 and 1.7, respectively (Woods et al. , 2000). These values are 



slightly smaller than the corresponding values for H Lyman a. We note that 
limb darkening or brightening can also affect the observed transit depths 



236 (KIO, Schlawin et al. , 2010). However, center-to-limb variations are not par- 



237 ticularly significant in the solar H Lyman ol and O I lines ( Curdt et al. 2008 



Roussel-Dupre , 1985) 



239 2.1.2. C II 1335 A and Si III 1206.5 A lines 



Linsky et al. (2010) used the medium resolution {R ~ 17,500) G130M 



241 grating of the COS instrument to observe four transits of HD209458b between 

242 September 19 and October 18, 2009 at wavelengths of 1140-1450 A. During 

243 each HST visit, they observed the star during transit, at first quadrature, 

244 secondary eclipse and second quadrature. In order to obtain the out-of- 

245 transit reference spectra, they co-added the secondary eclipse and quadrature 

246 observations from all four visits. They also co-added the in-transit spectra 
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from all visits to create a single spectrum. The detections of the transits in 
the C II 1335 A and Si III 1206.5 A lines were compared with observations 
of other lines such as Si IV 1395 A in which the transit was not detected. 



250 _We note that the transit in the Si III line was not detected earlier by Vidal 



Madjar et al. (2004). However, these authors use a wider wavelength interval 



to calculate the line-integrated transit depth and yet report a 2a upper limit 
of 5.9 %, implying that the 3(T upper limit includes the transit depth observed 



254 by Linsky et al. (2010). As we have seen, stellar variability can also cause 



changes in the perceived transit depth. Therefore we are not convinced of 
the reality of a non-detection in the earlier STIS G140L observations. 

In order to calculate transit depths in the C II and Si III lines, we created 



258 ijnodels for the stellar emission line profiles. We note that Ben-Jaffel and 



Hosseini (2010) also modeled the C II lines in order to fit the transit depth 



260 in the low resolution G140L data ( Vidal-Madjar et al. , 2004) where the main 



components of the C II multiplet are unresolved. Their model line profiles 
were also constrained by the high resolution STIS E140M observations of 



Vidal-Madjar et al. (2003) that resolve the components, although this data 



has low S/N and it was not used for transit observations. Here we use the 
higher S/N COS observations that also resolve the main components of the 
C II lines to constrain the line profile models. 

The atomic line parameters for the C II and Si HI lines are listed in 
Table [TJ The C II multiplet consists of three separate emission lines. The 
two lines at 1335.66 A and 1335.71 A (hereafter, the C II 1335.7 A line) are 
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1334.0 1334.2 1334.4 1334.6 1334.8 
Wavelength (A ) 



335.0 




1335.2 1335.4 1335.6 1335.8 1336.0 1336.2 
Wavelength (A ) 

Figure 1: (a) The C II 1334.5 A line of HD209458 (solid line, 
Linsky et al., 2010) fitted with a Voigt profile (see Table [l]) 
and adjusted for absorption by the ISM (dotted line). We as- 
sumed that the column density of ground state C+ in the ISM is 
2.23 X 10^^ m^^. The relative velocity of the ISM with respect 
to Earth is -6.6 km s^^ and the effective thermal velocity along 
the LOS to the star is 12.3 km s^^ (|Wood et al.) |2005D. (b) The 



C II 1335.7 A line of HD209458 fitted with a Voigt profile. Ab- 
sorption by the ISM was assumed to be negligible. The model 
profiles were convolved to a spectral resolution of R= 17,500. 



unresolved in the COS (and most solar) observations. The core of the C II 
1334.5 A line is strongly absorbed by the ISM whereas the C II 1335.7 A line 
is not similarly affected. The ground state ^Pi/2,3/2 of C+ is split into two 
fine structure levels. Interstellar absorption of the C II 1335.7 A emissions 
depends on the population of the ^^3/2 level in the ISM. In line with a similar 
assumption regarding the excited states of O (KIO), we assumed that this 
population is negligible. To estimate the emission line profile from HD209458, 



277 _we fitted the solar C II 1335.7 A line from the SUMER spectral atlas (Curdt 



et al. , 2001) with a Voigt function and adjusted the result to agree with 
observations (Linsky et al. , 2010). Figure [l] shows that the model line profile 
agrees reasonably well with the observed line profile. The fit parameters for 
this and other relevant emission lines are listed in Table [U 

Table 1: Atomic line parameters" 



Line 


Ao(A) 


fo 




s 






r 


F(10 erg cm ^ s ^) 


C II 


1334.53 


0.1278 


2.393 


X 


10^ 


0.12 


0.266 


1.16 


C IP 


1335.66 


0.01277 


4.773 


X 


10^ 


N/A 


N/A 


N/A 


C II 


1335.71 


0.1149 


2.864 


X 


10^ 


0.16 


0.224 


1.66 


Si III 


1206.5 


1.669 


2.550 


X 


10^ 


0.12 


0.535 


2.09 



Morton 



1991) 



"from 

*The line profiles are given by p{X) = [F/{AX£,^/^^)]V{a,u) where V is the IDL 
Voigt function with a = r/(47rAAD) and u = {X — Xo)/AXd- The flux is given at 
Earth distance (47 parsec). 

'^The emission line merges with the C II 1335.71 A line. 



282 Strong absorption by the ISM makes fitting the C II 1334.5 A line more 

283 ipomplicated than fitting the C II 1335.7 A line. Following Ben-Jaffel and 
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Hosseini (2010), we estimated that the column density of ground state C"*" 



in the ISM is 2.23 x 10^^ m ^ by scahng the column density of C"*" mea- 



286 sured along the LOS to Capella (Wood and Linsky, 1997) to the distance of 



HD209458. We fitted the solar C II 1334.5 A line from the SUMER spectral 
atlas with a Voigt profile, and used the estimated column density and the 



289 results of Wood et al. (2005) to calculate absorption by the ISM. We then 



290 varied the total flux within the line profile until the results agreed with the 



observations of Linsky et al. (2010). As a result, we obtained a pre-ISM flux 



ratio of [C II 1334.5 A]/[C II 1335.7 A] ~ 0.7, which agrees well with the solar 
value ( Curdt et al.[ 2001). The model and observed line proflles are again 
shown in Figure [Tj We note that the ISM is optically thick at wavelengths 
between 1334.54 A and 1335.58 A, which correspond to Doppler shifts of 2.2 
and 11.2 km s"^, respectively. However, the observed flux in this region is 
not zero because of spectral line broadening in the COS instrument. 

The C II lines are formed in the upper chromosphere and lower transi- 
tion region of the solar atmosphere. Similarly with H Lyman a, the brightest 



emissions are associated with plage activity (e.g., Athay and Dere, 1989 ). The 
plage and active network contrasts for the C II lines are 5.9 and 1.5, respec- 



302 tively (Woods et al. , 2000). The large spatial variability makes it difficult to 



characterize center-to-limb variations. However, different solar observations 
point to approximately 40 % limb brightening in the 1335.7 A line and prob- 
ably a similar variation in the 1334.5 A line, with the intensity rising steadily 



306 at (X > 0.6 (Lites et al. , 1978; Judge et al. , 2003). This brightening effect is 
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307 due to the broadening of the emission hne in the hmb. 

308 The Si III hne arises from a transition between the ground state and 

309 the excited state. Again, we fitted the Si III line profile from the SUMER 



310 spectral atlas (Curdt et al. , 2001) with a Voigt profile and adjusted the 



311 resulting line profile to agree with observations of IID209458 (Linsky et al. 



2010). Figure [2] shows the observed and model line profiles, and the fit 



313 parameters are listed in Table [Tj We assumed that the abundance of Si^"*" in 

314 the ISM is negligible. This is supported by a lack of detectable absorption 

315 by the ISM in the observed line profile. We note that absorption by the 

316 ISM affects the interpretation of the measured transit depths only if parts of 

317 the line profile are entirely absorbed or if the properties of the ISM change 

318 between observations. 

319 The solar Si III line has not been studied to the same degree as the C 

320 II lines. However, some constraints on the variability and center-to-limb 



321 variations of the Si III line have been obtained. For instance, Nicolas et al. 



322 (1977) used SKYLAB observations to study the center-to-limb variations of 

323 the Si III emissions from the quiet chromosphere. They found that the line 

324 is strongly limb-brightened. The total line intensity increased by a factor of 

325 2.4 between the disk center and /i = 0.73, and then again by a factor of 3 

326 towards the edge of the limb. This means that the total intensity in the limb 

327 is a factor of ~7 higher than at the disk center, and the line profile is also 

328 significantly broader at ;U > 0.73 with a self-reversal that does not appear at 

329 the disk center. The Si III emissions from the Sun also exhibit strong spatial 
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1206.0 1206.2 1206.4 1206.6 1206.8 
Wavelength (A ) 



1207.0 



Figure 2: 



Linsky et al 



The Si III 1206.5 A line of HD209458 (solid line, 
2010) fitted with a Voigt profile (dotted hne, see 



Table T]). Absorption by the ISM was assumed to be neghgi- 
ble. The model profile was convolved to a spectral resolution of 
R = 17,500. 
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330 and temporal variability (e.g., Nicolas et al. 1982). Limb brightening makes 



331 the transit depth appear smaller when the planet is covering the stellar disk 



332 while a steeper transit is seen during ingress and egress (e.g., Schlawin et al. 



333 2010). On the other hand, if the planet covers active regions on the disk, the 

334 transit depth can appear significantly deeper (KIO). 

335 Spatial variability and limb brightening of the emission lines on HD209458 

336 can potentially be studied through a careful analysis of the transit light 

337 curves. Ideally, the observations should be analyzed in the time tag mode 



338 (e.g., Ben-Jaffel, 2007) to identify variations. This type of reanalysis of the 

339 COS data is beyond the scope of this paper. Instead, we use idealized transit 

340 depths based on a uniformly emitting stellar disk in Section [3] to show that 

341 the optical depth of the extended thermosphere in the FUV lines is significant 

342 and that the transit in the O I, C II, and Si III lines is comparable to the 

343 transit in the H I line. We consider this sufficient for the present purposes. 

344 2.2. Empirical model 

345 KIO developed an empirical model to fit the UV transit observations of 

346 HD209458b and other extrasolar planets. They argued that the H I tran- 

347 sits can be explained in terms of three simple parameters that describe the 

348 distribution of H in the thermosphere. These parameters are the pressure 

349 Pq where II2 dissociates (the bottom of the H layer), the mean temperature 

350 T within the H layer and an upper cutoff level based on the ionization 

351 of H. Transits in other emission lines can be explained by fitting the abun- 
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352 dance of the heavier absorbers with respect to H. The details of this model 

353 are discussed by KIO and are not repeated here. Basically it calculates the 

354 transit depth observed at Earth orbit by assuming that the planet and its 

355 atmosphere constitute a spherically symmetric obstacle with a density profile 

356 in hydrostatic equilibrium up to the sonic point (when a sonic point exists). 

357 Absorption by the ISM and spectral line broadening within the observing 

358 instrument (STIS or COS) are taken into account. In Section [3| we compare 

359 the transit depths calculated by the empirical model with results based on 

360 the hydrodynamic model (see below) to show how the empirical model can 

361 be used to fit the data and guide the development of more complex models 

362 of the upper atmosphere. 

363 2.3. Hydrodynamic model 

364 We developed a one- dimensional, non-hydrostatic escape model for HD209458b 

365 to constrain the mean temperature and ionization in the upper atmosphere 

366 (Paper I). Results from this model demonstrate that the empirical model 

367 is physically meaningful. It solves the vertical equations of motion for an 

368 escaping atmosphere containing H, H+, He, He"'', C, C+, O, O"*", N, N+, 

369 Si, Si+, Si^"*", and electrons. The model includes photoionization, thermal 

370 ionization, and charge exchange between atoms and ions. It calculates the 

371 temperature profile based on the average solar X-ray and EUV (XUV) flux 



372 jp,nd up to date estimates of the photoelectron heating efficiencies (Cecchi 



373 Pestellini et al. , 2009, Paper I). The lower boundary of the model is at 1 fihar 
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374 and thus the model does not include molecular chemistry. This is justified 

375 because molecules are dissociated by photochemical reaction networks near 



376 the 1 /ibar level (e.g., Garcia Munoz 2007 Moses et al. 2011). The upper 



377 boundary of the model is typically at 16 Rp. We placed the upper boundary 

378 at a relatively high altitude above the region of interest in this study, which 

379 is below 5 Rp. However, we do not consider the results to be necessarily 

380 accurate above 5 Rp (see Paper I for further details). 

381 3. Results 

382 3.1. Transit depths 

383 In this section we constrain the temperature and composition of the upper 

384 atmosphere of HD209458b through a combined analysis of transit observa- 

385 tions in the FUV emission lines. We also compare results from a hydrody- 



386 namic model (see Section 2.3 and Paper I) with the empirical model of KIO, 

387 and confirm that the latter can be used to constrain the basic properties of 

388 the density profiles in the thermosphere. 

389 3.1.1. Neutral atoms 

390 KIO demonstrated that absorption by hydrogen in the extended thermo- 

391 sphere of HD209458b explains the transits in the H Lyman a line, and used 

392 the observations to constrain the mean temperature and composition of the 

393 thermosphere. They fitted both the line-integrated transit depth based on 



394 iJ^he low resolution G140L data ( Vidal-Madjar et al. 2004 Ben-Jaffel and 
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Hosseini, 2010), and the transit depths and hght curve based on the medium 



2003; 


Ben- Jaffel 


2007 


2008 



400 r, 



(see Figures 5 and 6 of KIO for the results). The results imply that the lower 
boundary of the absorbing layer of H is at Pq = 0.1-1 /ibar, the mean temper- 
ature within the layer is T = 8,000-11,000 K, and the upper boundary is at 
= 2.7 Rp. Recent photochemical calculations imply that H2 dissociates 



near the 1 /ibar level [e.g.. Paper I, Moses et al. (2011)], and this is also 
supported by an observational lower limit for the vertical column density of 



H (France et al. , 2010). Hence the mean temperature in the thermosphere 
of HD209458b is approximately 8,250 K (the M7 model of KIO). 

Here we show from simple arguments that the M7 model agrees in princi- 
ple with the observed H Lyman a transit depth. A similar procedure can also 
be used to obtain crude estimates of the transit depth for other systems, pro- 
vided that the properties of the line profile and the ISM are known. Later in 
this section we compare the results of the empirical model with results from 
the hydrodynamic model, and show that they are consistent. To start with. 



iJ;he line-integrated H Lyman a transit depth of 6.6 ± 2.3 % (Ben- Jaffel and 



Hosseini 



2010) is consistent with a 6 % transit depth at 1215.2 A (a Doppler 
shift of -120 km s~^ from the line core) i.e., the blue peak of the observed 
line profile (Figure 5 of KIO). Assuming that the extended atmosphere is 
spherically symmetric, a 6 % transit depth measured in the blue wing (bw) 
of the line profile implies an optical depth of Tbw ~ 1 at the impact parameter 
of 2.1 i?p from the center of the planet. Figure [3] shows the absorption cross 
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section of H in the Lyman a line at a temperature of 8,250 K. The cross 
section at 1215.2 A is dbw = 2 x 10^^^ m^ and thus a LOS column density of 
A^H = 5 X 10^^ m^^ at 2.1 Rp is required to explain the observed absorption. 



Ekenback et al. (2010) and Lammer et al. (2011) argued that the optical 



depth of H in the thermosphere of HD209458b is not significant and instead 
a large cloud of energetic neutral atoms (ENAs) is required to explain the 



H Lyman a observations. In particular, Lammer et al. (2011) claimed that 
a column density of A^h ~ 10'^^ m~^ in the thermosphere is required for 
strongly visible absorption. It is easy to see that this estimate is not correct 
because the wings of the line profile become optically thick with column 



428 densities much smaller than this. It also disagrees with Garcia Munoz (2007) 



429 and Ben-Jaffel (2007, 2008) who were the first to suggest that H in the 



thermosphere may be sufficiently abundant to explain the observations. This 



431 |j3asic result has also been confirmed by more recent calculations by Trammell 



et al. (2011). All of these calculations show that the optical depth of the 



thermosphere below 3 Rp is not negligible. 

The empirical model of KIO is simplified by the use of the hydrostatic 
approximation. This is well justified even if the atmosphere is escaping. In 
general, the density profile of the escaping gas in the thermosphere can be 



437 estimated from (Parker, 1964): 



exp 



^ du dv 
I du 



(1) 
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Figure 3: The absorption cross section of H in H Lyman a 
line at T = 8,250 K. Thermal and natural line broadening were 
modeled with a Voigt profile. 

438 where ^ = r/ro, c^(0 = kT{$^)/m, W = GMp/ro, m is the molecular weight, 

439 and V is the vertical flow speed. For convenience, we retained Parker's original 

440 notation. The first integral on the right hand side of equation ([T]) applies in 

441 hydrostatic equilibrium. The second integral is negligible below the sonic 

442 point for transonic escape and always negligible for subsonic escape (i.e., 

443 evaporation). The sonic point on HD209458b is always above 3 Rp (Paper I) 

444 and thus the density profile of H is approximately hydrostatic at least up to 

445 2.1 Rp. 

446 Assuming a hydrostatic atmosphere, the LOS column density of H at 2.1 
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447 Rp can be estimated from: 



Nii{r = 2.1 Rp) ^ rioexp 



GMpm (\ 



kT 



r To 



(2) 



448 where T and H are the mean temperature and scale height, respectively. We 

449 note that the mean thermal escape parameter is: 



X(r) = 

kTr 



(3) 



Assuming that T = 8,250 K and using the planetary parameters of HD209458b 
(Mp = 0.7 Mj, Rp = 1.3 Rj), we obtain values of Xq = 13.8 and X(2.1 Rp) = 
6.6. Thus, according to equation the column density of = 5 xlO^^ 
m~^ that is required to explain the observations implies that no = 3.5 x 10^^ 
m^'^. This in turn means that po = 0.4 yubar i.e., close to the lower boundary 
of the M7 model (here the agreement with the M7 model is obviously not 
exact because the transit depths in KIO are based on a complete forward 
model of the observed transit within the whole line profile). 

The parameters of the empirical model can be compared with correspond- 
ing values derived from the hydrodynamic models presented in Paper I (see 



Section 2.3 here for a brief summary). For instance, the mean temperature 
of the empirical model corresponds roughly to the pressure-averaged temper- 



462 ature of the thermosphere, which is given by (e.g., Holton, 2004): 



T 



/;;T(p) d(lnp) 



\n{p2/pi 



(4) 
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463 The hydrodynamic calculations show that the pressure averaged (mean) tem- 

464 perature below 3 Rp based on the average solar flux varies between 6,000 K 

465 and 8,000 K for net heating efficiencies rj^et between 0.1 and 1. This tempera- 

466 ture is relatively insensitive to different assumptions about heating efficiencies 

467 or the upper boundary conditions. In the reference C2 model of Paper I the 

468 mean temperature is 7,200 K. This model is based on our best estimate of 

469 the heating efficiencies that are appropriate in the strongly ionized upper 

470 atmosphere of HD209458b. 

471 With a cutoff level at 2.7 Rp, we obtained a line-integrated H Lyman a 

472 transit depth of 4.7 % based on the density of H in the C2 model. This 



value agrees with the observations to within la ( Vidal-Madjar et al. , 2004 



Ben-Jaffel and Hosseini, 2010), but it is smaller than the transit depth of 



475 6.6 % predicted by the M7 model. One reason for this is the lower mean 

476 temperature of the C2 model. In order to facilitate a direct comparison 

477 between the hydrodynamic model and the empirical model, we calculated the 

478 empirical transit depth based on the mean temperature of 7,200 K (hereafter, 

479 the M7b model). The line-integrated transit depth based on this model is 

480 5.8 %, which is still higher than the transit depth based on the C2 model. 

481 Figure |4] shows the density profiles of H, H"*", O, and O"^ from the C2 

482 model, and the density profile of H from the M7b model. The difference 

483 between the transit depths based on the empirical and hydrodynamic models 

484 arises because the C2 model has large temperature gradients (Paper I) and a 

485 gradual H/H"*" transition rather than a sharp cutoff. The difference does not 
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12 5 4 5 

Figure 4: Density profiles O and H based on the C2 model 
(Paper I) and the density profile of H based on the empirical 
model of KIO with a mean temperature of 7,200 K. 

486 arise because the density profile of tlie C2 model deviates from hydrostatic 

487 equilibrium. Given the temperature gradient in the model, the density profile 

488 is almost exactly in hydrostatic equilibrium below 3 Rp. In fact, the neutral 

489 density profile of the C2 model is better represented by a mean temperature 

490 of 6,300 K (not shown). This implies that the correspondence of the mean 

491 temperature of the empirical model and the pressure averaged temperature 

492 of the hydrodynamic model is relatively good but not exact. 

493 A better agreement between the transit depths based on the C2 and M7b 

494 models is obtained with the cutoff level of the C2 model at 5 Rp (see Table [2|. 

495 Figure [5] shows the line-integrated transit depth within the H Lyman a line 

496 profile as a function of the cutoff level for the C2 model. This figure indicates 
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Figure 5: Line-integrated transit depth within the H Lyman a 
line as a function of the obstacle cutoff level based on the C2 
hydrodynamic model (Paper I). The KIO density profile applies 
to the M7b model (see text). 



497 that the transit depth increases less steeply with altitude above 4 Rp than 

498 it does below this cutoff level, and saturates near 5 Rp. This is a natural 

499 consequence of the fact that the LOS column density decreases approximately 

500 exponentially with altitude in the lower thermosphere and the wings of the 

501 line profile become optically thin at high altitudes. 

502 We note that KIO already compared the results based on the M7 model 

503 with in-transit transmission as a function of wavelength within the H Ly- 



504 ipan a line profile and the light curve derived from the G140M data (Ben 



Jaffel, 2007, 2008). It is not necessary to repeat a similar comparison here 



because spherically symmetric models that predict line-integrated transit 
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depths that agree with the measured values are generally compatible with 
both the G140L and G140M data. This is partly because the uncertainty of 
the individual data points within the line profile and the light curve is large, 
and thus they do not strongly constrain the properties of the atmosphere. 

According to Figure |4| the H/H"*" transition in the C2 model occurs near 
3.1 Rp. The exact location depends on photochemistry and vertical velocity, 
and generally the transition occurs near or above 3 Rp (Paper I). With a 
fixed pressure at the lower boundary, a faster velocity leads to a transition 



a higher altitude. These results disagree with Yelle (2004) and Murray- 



Clay et al. (2009) who predicted a lower transition altitude, but they agree 



517 qualitatively with the solar composition model of Garcia Munoz (2007). The 



density profiles of O and are strongly coupled to H and H"*" by charge 
exchange. As a result, the 0/0+ transition occurs generally near the H/H"*" 
transition. For instance, in the C2 model it is located near 3.4 Rp. We note 
that significant ionization of H and O above 3 Rp is anticipated by KIO and 
there is good agreement on this between the empirical and hydrodynamic 
models. 

The detection of O at high altitudes constrains the mass loss rate and 



the ionization state of the upper atmosphere (see Section 3.2). However, 
the large uncertainty in the observations means that repeated observations 
are required to confirm the transit depth. The M7 model of KIO with a 



solar 0/H ratio (Lodders, 2003) yields an O I transit depth of 3.9 %, which 
is within 1.5a of the observed value and therefore a satisfactory fit to the 
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530 data. This value agrees well with the O I transit depth based on the C2 

631 model (see Table |2]). KIO argued that a higher transit depth is possible if the 

532 mean temperature is higher and/or if the 0/H ratio is enhanced with respect 

533 to solar. The hydrodynamic calculations indicate that the latter option is 

534 more favorable because higher temperatures lead to stronger ionization of 

535 O and may not help to significantly enhance the transit depth. Indeed, the 

536 predicted transit depth agrees with the observations to better than la if the 

537 0/H ratio is enhanced by a factor of 5 relative to solar (see the MS0L2 

538 model in Table I2I). 



Model 



Table 2: Model parameters and transit depths (%) 



M (10' kg s--') rjnet" 



(K) 



H I"^ 
6.6 ± 2.3=' 



01 C II 1334.5 A C II 1335.7 A 

10.5 ± 4.4 7.6 ± 2.2 7.9 ± 1.5 



Si III 
1.2 ± 1.4 



CI 
C2 
C3 
C4 
C5 
C6 
SOL2 
MSOL2 
M7 



5.6 
4.0 
6.4 
4.5 
5.6 
3.9 
11.0 
6.0 
N/A 



0.56 
0.44 
0.57 
0.46 
0.56 
0.45 
0.50 
0.66 
N/A 



7250 
7200 
6450 
7110 
7290 
7310 
7390 
7370 
8250 



5. 5 
5.7 
5.4 
5.2 
5.6 
5.9 
5.0 
5.6 
6.6 



4.1 
4.0 
3.7 
3.7 
4.2 
4.1 
4.2 
7.1 
3.9 



3.4 
3.2 
3.1 
3.2 
3.5 
3.3 
4.3 
5.5 
3.9 



6.8 
6.7 
5.5 
6.3 
6.9 
6.9 

r.8 

10.5 
8.0 



5.0 
4.6 
3.5 
4.6 
5.1 
4.6 
6.8 
8.0 
5.8 



"Net heating efficiency (see Section 2.3) i.e., the ratio of the net heating flux at all 
wavelengths to the unattenuated stellar flux (0.45 W m~^) at wavelengths shorter 
than 912 A. 

^Pressure averaged temperature below 3 Rp. 

"^Line-integrated transit depth. The upper boundary of the absorbing atmosphere 
is at 5 Rp apart from the M7 model where it is at 2.7 Rp for neutrals and at 5 Rp 
for ions. 

"^Observed values from 



Ben-Jaffel and Hosseini 



(2010) and Linsky et al. (2010) 



with la errors. 



539 We have now verified that the empirical model can be used to constrain 

540 the mean temperature and extent of the absorbing layer in the thermosphere 
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541 of HD209458b. In particular, the comparison of the empirical model with 

542 the hydrodynamic model shows that the results of KIO were not affected 

543 by the simplifying assumption of hydrostatic equilibrium. We note that the 

544 purpose of the empirical model is to identify physical processes that might 

545 otherwise be missed in more complex models that are often based on a large 

546 number of uncertain assumptions. The results from any model can now be 

547 compared with the observations by identifying the limits of the absorbing 

548 layer and calculating the global pressure averaged temperature within that 

549 layer. The values can then be compared with the parameters of the best-fit 

550 empirical model. 

551 3.1.2. Ions 

552 Figure [6] shows the density profiles of the carbon and silicon ions in the 

553 thermosphere of HD209458b based on the C2 hydrodynamic model (Paper 

554 I). The major carbon and oxygen-bearing species in the lower atmosphere, 

or 



555 CO and H2O, dissociate near the 1 /ibar level (e.g., Moses et al. , 2011 



556 Lavvas et al., in preparation). Thermochemical calculations indicate that 



557 SiO is the dominant silicon-bearing gas on HD209458b ( Visscher et al. , 2010 ). 

558 The detection of Si^"*" in the upper atmosphere implies that the formation of 



559 silicon clouds in the lower atmosphere is suppressed (see Section 3.3), and 

560 SiO is also dissociated near 1 /ibar either thermally or by photochemistry. 

561 Thus we assumed that only atomic carbon and silicon are present in the 



562 thermosphere, initially with solar abundances (Lodders, 2003). 
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Figure 6: Density profiles of C and Si in the upper atmosphere 
of HD209458b based on the C2 model (Paper I). 

563 According to Figure [6| the C/C^ transition occurs at a much lower ahi- 

564 tude of 1.2 Rp than the H/H+ and 0/0+ transitions. Sihcon is also almost 

565 fully ionized at the lower boundary of the model, and the Si^"'"/Si"*" ratio is 

566 about 10 % below 3 Rp. These results are in qualitative agreement with the 

567 observations, because they show that H and O are mostly neutral below 3 Rp 

568 whereas C and Si are mostly ionized, and a significant abundance of Si^"*" is 

569 possible. However, it is also useful to explore if the results from the models 

570 are in quantitative agreement with the observations and if not, to adjust the 

571 model parameters as necessary to explain the data. 

572 With a cutoff level at 2.7 -Rp, we obtained line-integrated transit depths 

573 of 2.3 %, 3.6 %, and 2.5 % in the C II 1334.5 A, C II 1335.7 A and Si III lines, 

574 respectively, based on the C2 model. Here we assumed that the population 
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675 of the levels of C"*" are in LTE. We also calculated empirical transit depths 

676 based on the M7 model. In order to do this we assumed that both C and 

677 Si are ionized at the lower boundary, and that 10 % of silicon is Si^'''. With 

678 these assumptions we obtained C II 1334.5 A and 1335.7 A transit depths of 

579 2.7 % and 4.2 %, respectively, and a Si III transit depth of 3 %. These values 

580 agree well with the C2 model, and further demonstrate the consistency of 
681 the empirical and hydrodynamic models. However, they deviate from the 

582 observed values by more than 2a. 

583 The cutoff level of the empirical model is somewhat arbitrary. For neutral 

684 species it is partly based on ionization (KIO), but this criterion obviously 

685 does not apply to ions. With a cutoff level at 5 Rp for the ions only, the 
586 M7 model yields line-integrated transit depths of 3.9 %, 8 %, and 5.8 % in 

687 the C II 1334.5 A, C II 1335.7 A, and Si III lines, respectively, if 40 % of 

688 silicon is Si^"*". These values agree with the observed values to better than 
589 2a. Similarly, by extending the cutoff level of the C2 model to 5 Rp, we 

690 obtained transit depths of 3.2 %, 6.7 %, and 4.6 % in the C II 1334.5 A, 

691 C II 1335.7 A, and Si III lines, respectively (see Table [2|. These values 
592 deviate from the observed values by 2a, 0.9a, and 2. 60", respectively. The 

693 transit depths predicted by the M7 model are higher partly because the mean 

694 temperature of 8,250 K is higher than the corresponding temperature in the 
595 C2 model (Paper I). This also leads to the higher Si^^/Si+ ratio that we used 

696 here. 

697 It is not clear if the apparent disagreement between the models and some 
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598 of the observations needs to be taken seriously. Stellar activity and other un- 
699 certainties mean that the true transit depth can differ from measured values 



600 by a significant factor (see Section 2.1). Further, the C2 model agrees with 

601 the line-integrated H Lyman a and C 11 1335.7 A transit depths to within la, 

602 and with the O 1 and C 11 1334.5 A lines to within 2o. Thus we could argue 

603 that the present observations are roughly consistent with solar abundances 

604 and heating based on the average solar XUV flux. Nevertheless, we explore 

605 the apparent disagreement between the C2 model and the observations fur- 

606 ther below. This disagreement is limited to the O 1, C 11 1334.5 A, and Si 

607 111 lines. 

608 Figure [7] shows the observed in-transit flux differences in the C 11 and Si 

609 111 lines as a function of wavelength together with different model predictions. 

610 The observations indicate that the transit depths based on the C2 model fall 

611 short of the observed values because the model underestimates the width 



612 of the absorption lines. Linsky et al. (2010) argued that there is velocity 

613 structure within the line proflles near Doppler shifts of -10 km s~^ and 15 

614 km s~^ that accounts for broad absorption. However, the uncertainty of the 

615 individual data points is too large to constrain the shape of the absorption 



616 lines in detail. We agree with Linsky et al. (2010) that the presence of velocity 

617 structure needs to be conflrmed by future observations. Thus the additional 

618 absorption could also arise from spectral line broadening. 

619 We note that H2 absorbs within the Si 111 line and it can contribute 

620 to the observed absorption. The cross section has a lot of structure in this 
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Figure 7: Flux differences between tfie stellar in-transit and 
out-of-transit C II and Si III emission lines. The data points 



were taken from Linsky et aU^2010). Also shown are model 



predictions based on the C2 model (dashed hne) and the S0L2 
model (dash-dotted line). See text and Table [2] for the details 
of these models. 



wavelength region, but for the sake of the argument we adopted a mean cross 



622 section of ~2 x 10^^'^ m^^ (Backx et al. , 1976). A transit depth of 8 % across 

623 the hne profile implies that the atmosphere is optically thick up to 2.4 Rp 

624 where the LOS optical thickness is r ~ 1, and thus a H2 column density of 

625 5 X 10^^ m~^ along this LOS would be required to explain the observed transit 

626 depth. This is unrealistic because the required column density is higher than 



627 the corresponding column density of H (see Section 3.1.1). Also, according 

628 to our photochemical calculations H2 dissociates in the thermosphere, and 

629 the mixing ratio of H2 falls below 0.1 above the 0.1 /ibar level. Such a low 

630 abundance of H2 has no effect on the Si III transit depth. 

631 A higher mean temperature leads to higher transit depths. Therefore we 

632 generated a new hydrodynamic model by multiplying the average solar flux 

633 by a factor of 2 and assumed a net heating efficiency of r^net = 0.5 (here- 

634 after, the S0L2 model). This model agrees to better than la with all of the 

635 observed transit depths apart from the O I and C II 1334.5 A lines (see Fig- 

636 ure[7]and Table [2]). The mean temperature of the model below 3 -Rp is 7,400 

637 K, which is lower than the mean temperature in the M7 model. However, 

638 absorption by the S0L2 model is strengthened by velocity dispersion within 

639 the escaping plasma that is not included in the M7 model. In general, the 

640 outflow velocity of the S0L2 model is significantly higher than the velocity 

641 in the C2 model. The model also predicts a mass loss rate of 10® kg s~^, 

642 which is twice as high as the mass loss rate based on the C2 model. This 

643 proves that a higher stellar XUV fiux or a corresponding alternative energy 
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644 source can explain the observations. An enhancement of the average solar 

645 flux by a factor of 2 is not unreasonable, and would roughly correspond to 

646 solar maximum conditions. 

647 In Paper I we noted that the stellar XUV flux, or the corresponding 

648 alternative heat source, would have to be 5-10 stronger than the average 

649 solar flux to produce a mean temperature between 8,000-9,000 K. Under 

650 such circumstances, the predicted transit depths in the C II and Si III lines 

651 would obviously be even higher than the values predicted by the S0L2 model. 

652 Indeed, higher temperatures broaden absorption in the wings of the line 

653 proflles and may help to explain the in-transit flux differences better (see 

654 Figure [T]). However, the energy input and temperature in the model cannot 

655 be increased without bound. Higher temperatures and flux lead to more 

656 efficient ionization of the neutral species, and as a result the transit depths 

657 in the H Lyman a and O I lines begin to decrease. Also, mass loss rates of 

658 1 0^-10^° kg s~^ lead to the loss of 10-100 % of the planet's mass over the 

659 estimated lifetime of the system, and this probably limits reasonable energy 

660 inputs to less than ~10 times the solar average on HD209458b. 

661 In addition to higher temperature and velocity, supersolar abundances of 

662 O, C, and Si can also lead to higher transit depths. This option is interesting 

663 because it also allows for a higher transit depth in the O I lines. As an 

664 example, we generated the MS0L2 model by enhancing the solar 0/H, C/H, 

665 and Si/H abundances in the hydrodynamic model by a factor of 5. As a 

666 result, we obtained transit depths that agree with nearly all of the observed 
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line- integrated transit depths to better than la (see Table [2]). However, the 
MS0L2 model overestimates the line-integrated C II 1335.7 A transit depth, 
and generally overestimates absorption within the cores of the C II and Si 
III lines. This could imply that a higher temperature or some other source 
of additional broadening is a better explanation of the Si III and C II transit 
depths while a supersolar O/H ratio is required to match the measured O I 
transit depth. However, the data points in Figure [7] and the observed O I 
transit depth are too uncertain and do not allow for firm constraints on this. 

Enrichment of heavy elements is a common feature on the gas and ice 
giants in the solar system. For instance, the C/H, N/H, S/H, Ar/H, Kr/H, 
and Xe/H ratios in the atmosphere of Jupiter are all enriched by factors of 



678 2-3 with respect to solar abundances (e.g., Mahaffy et al. , 2000; Wong et al. 



2004). On Saturn, on the other hand, the C/H ratio is enriched by a factor 



of 10 (Flasar et al. , 2005 Fletcher et al. , 2009). Enrichment by factors of 



4-20 is expected in the N/H and S/H ratios, although condensation of NH3 
and H2S in the deep atmosphere of Saturn makes it difficult to constrain the 



abundances precisely (see Fouchet et al. , 2009, for a review). On Neptune 
and Uranus the C/H ratio is believed to be enriched by factors of 30-50 



685 (e.g., Owen and Encrenaz 2003 Guillot and Gautier 2007), and similar 



enrichments are possible in the abundance ratios of some of the other heavy 
elements. Substantial enrichment of heavy elements with respect to solar 
abundances is therefore also feasible in EGP systems even if the metallicity 
of the star is close to solar. 
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690 Unfortunately, we cannot use the current observations to constrain the 

691 elemental abundances of the atmosphere with accuracy. In this regard, the 

692 large uncertainty of the observations is unfortunate, because similar observa- 

693 tions can potentially be used to estimate them. The dissociation of molecules 

694 at the relatively high pressure of 1 /ibar and the lack of diffusive separation 

695 mean that the abundances of the heavy atoms and ions are simply depen- 

696 dent on the elemental abundances and ionization rates. Observations of the 

697 neutral species can therefore be used to constrain the temperature and ion- 

698 ization state, and thus the elemental abundances of the heavy species, but 

699 the S/N of the current data does not allow for strong constraints. 

700 It is interesting to note that while the velocity structure of the escap- 

701 ing plasma can lead to broader absorption that helps to explain the transit 

702 depths, it is not necessarily detectable in the data. For instance. Figure [7] 

703 shows the transit depths based on the S0L2 model that has a relatively high 

704 radial velocity reaching 11 km s^^ by 5 Rp. The velocity structure is not 

705 detectable because the optical depth of the high velocity material is not suf- 

706 ficient, the LOS velocity at the limb of the planet is slower than the radial 

707 velocities in general, and because spectral line broadening within the COS 

708 instrument smooths the structure out of the line profiles. If the presence of 



709 velocity structure is confirmed in the data (Linsky et al. , 2010), it probably 

710 implies that there is detached, optically thick plasma moving at large veloc- 

711 ities around the planet. If this turns out to be the case, interaction with the 

712 stellar wind probably plays a role in giving rise to the observed absorption. 
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713 Such interaction may also produce turbulence that can broaden the absorp- 



714 tion further (e.g., Tian et al. , 2005). However, we note that non-thermal 



715 broadening such as that proposed by Ben-Jaffel and Hosseini (2010) does 



716 not appear to be necessary to explain the current observations. 

717 3.2. Ionospheric escape 

718 The escape of heavy atoms and ions has interesting consequences for the 

719 nature of the upper atmosphere. Here we discuss these consequences based on 

720 simple analytic arguments, and without explicit use of any complex models. 

721 The detection of heavy neutral species can be used to constrain the mass 

722 loss rate while the detection of heavy ions outside the atmosphere of the 

723 planet potentially constrains the strength of the planetary magnetic field. 



724 For instance, Hunten et al. (1987) derived an expression for the crossover 

725 mass limit for a neutral species s to be dragged along by an escaping 

726 neutral species t of mass nit < ms'. 

Tnc = mt^ 2 (5) 

727 where Ft is the flux (s~^ sr~^) of species t, xt is the volume mixing ratio, go is 

728 gravity at the lower boundary of the model region, and the mutual diffusion 

729 coefficient can be roughly estimated from: 

nDst = 1.52 X 10^^ \ ^ + ^] Vf cm-^s"^ (6) 
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730 where the masses M are in amu. We used equation ([s]) to estimate the mass 

731 loss rate that is required to drag neutral O to the exosphere of HD209458b. 

732 Assuming that Xt ~ 1, (^o = 10 m s~^, tq = Rp, T = 7,200 K, and nDst = 

733 1.3 X 10^^ m~^s~^, we obtain Ft ^ 2.8 x 10^^ s~^ sr~^. This implies a 

734 minimum mass loss rate of 6 x 10^ kg s~^. 

735 The ionosphere of HD209458b is mostly neutral below 3 Rp but even weak 

736 ionization can lead to frequent Coulomb or ion-neutral collisions that enable 

737 heavy ions or atoms to escape more efficiently. In order to illustrate the role 

738 of different collisions in transporting O and Si"*", Figure [s] shows the collision 

739 frequencies for these species with H and H"^ as a function of altitude based 

740 on the C2 model (Paper I). We used approximate expressions for resonant 

741 i^nd non-resonant ion- neutral collisions, and Coulomb collisions from [Schunk 
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and Nagy (2000) to calculate the momentum transfer collision frequencies. 

743 The collision frequency between two neutral species, on the other hand, was 

744 estimated from the mutual diffusion coefficient as: 



1 1 y'^^ntVT 



= 5.47 X 10-- ( - + - ) ^ (7) 



745 where the number density rtt is in cm~^ . The results indicate that the trans- 

746 port of O depends on collisions with H below 3.5 -Rp whereas the transport 

747 of Si+ depends on collisions with at all altitudes. 

748 Oxygen is the heaviest neutral species detected in the escaping atmo- 

749 sphere, and this implies that the mass loss rate from HD209458b is M > 6 x 
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Figure 8: Momentum transfer collision frequencies based on the 
C2 model. 



10^ kg s ^. This result agrees with Vidal-Madjar et al. (2003) although it 



is less model-dependent and based on different criteria. The dominance of 
Coulomb collisions means that the heavy ions can escape even if the mass loss 
rate is lower than this. Our models predict mass loss rates of M 5 x 10^ 
kg s^^ (Table [2]) and thus diffusive separation does not take place in the ther- 
mosphere of HD209458b for neutral species with masses less than ~130 amu. 
We note that this is the case even if escape is subsonic. In fact, equation ([s]) 
was originally derived for subsonic escape under the diffusion approximation 



although it is also valid for supersonic escape (Zahnle and Kasting, 1986). 

We used the collision frequencies to derive expressions for the ion fractions 
fi = at which ion-neutral and Coulomb collisions become important. 

The ratio of the non-resonant neutral-ion to neutral-neutral collisions exceeds 
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762 10 % when 



763 where i denotes H+, s denotes the colhding species, 7^ is the neutral po- 

764 larizabihty and all units are in cgs. However, the collision of O with H''" is 

765 resonant and in this case the required ratio differs slightly from the above 

766 expression. The ratio of the Coulomb to non-resonant ion-neutral collisions, 

767 on the other hand, exceeds 10 % when 



T3/2 M, 



768 where i denotes (or the dominant ion) and n is H (or the dominant 

769 neutral). For Si"*" this fraction is fi ~ 10~^ (with 7h = 6.7 x 10"^^ cm'^). 

770 These equations can be used to determine if equation ([s]) is valid, or if more 

771 complex plasma models are required. 



772 Trammell et al. (2011 ) argued that HD209458b could have a strong plane- 

773 tary magnetic field that can impede the escape of ions from equatorial regions 

774 and restrict it to the polar regions. Although the magnetic field does not di- 

775 rectly interfere with the escape of the neutral atoms, the trapped ions can 

776 stop them from escaping if the neutral-ion collision frequency is sufficiently 

777 high. Unfortunately, transit observations are not spatially resolved and they 

778 cannot be used directly to determine if escape is limited to the poles, or if the 

779 atmosphere is also escaping over the equator. However, the transit depths 

780 depend on the size of the optically thick obstacle covering the star. If mass 
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781 loss is suppressed at low and mid-latitudes, the heavy species are no longer 

782 mixed into the upper atmosphere other than at the poles where they are 

783 allowed to escape. This means that the cross-sectional area covered by the 

784 ions shrinks and may become insufficient to explain the observations. Even if 

785 the plasma spreads to cover a larger area after being ejected from the poles, 

786 it is diluted in the process and thus it is not clear if the resulting cloud would 

787 have sufficient optical depth to be detectable. 

788 Assuming that the charged particles escape the Roche lobe of the planet 

789 unimpeded at the equator, we estimated an upper limit for the planetary 

790 magnetic field by evaluating the magnetic moment that allows them to do 

791 so. In order to estimate this limit, we calculated the plasma /3 and the inverse 

792 of the second Cowling number (Co^^) below 5 Rp in the C2 model from: 

793 where v is the vertical velocity. Assuming a dipolar magnetic field, we ob- 

794 tained a limiting magnetic moment of 3.2 x 10^^ Am^ or 0.04 mj {mj = 

795 1.5 xlO^'' Am^ is the magnetic moment of Jupiter). This magnetic moment 

796 ensures that /3 > 10 below 5 Rp and that Co~^ reaches 10 by 5 Rp. Magnetic 

797 moments of rrip < 0.04 mj agree quite well with the scaling laws discussed 



798 by Griesmeier et al. (2004). We note that the limiting moment produces an 

799 equatorial surface field that is approximately 4 times lower than the surface 

800 field of the Earth. 
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3.3. Cloud formation on HD 20 945 8b 

We have shown that a substantial abundance of sihcon in the upper at- 
mosphere is required to produce a detectable transit in the Si III line. If the 
silicon ions originate from the atmosphere of the planet, at least a solar Si/H 
ratio is necessary. According to thermochemical equilibrium models, silicon 
should condense into clouds of forsterite (Mg2Si04) and enstatite (MgSiOa) 



in the lower atmosphere of HD209458b (Visscher et al. , 2010). If the forma- 



tion of enstatite is suppressed, silicon should condense to form quartz (Si02) 
instead. In any case, condensation is expected to remove almost all of the 
silicon from the upper atmosphere. The detection of Si^"*" implies that the 
abundance of silicon in the thermosphere is substantial, and thus the conden- 
sation of silicon does not take place in the atmosphere of HD209458b. This 
has significant implications for the structure and dynamics of the atmosphere. 
Sing et al. ( 2008a|[b ) analyzed the absorption line profile of Na in the 



atmosphere of HD209458b in detail and argued that the abundance of Na 
is depleted above the 1 mbar level. They suggested that this is due to the 
condensation of sodium into Na2S, although ionization could not be ruled 
out decisively. Based on the condensation temperature of Na, they argued 
that the temperature in the upper atmosphere of HD209458b near 1 mbar 
is 420 ± 190 K. This temperature is significantly lower than the outcome of 



typical radiative transfer models for HD209458b (e.g.. Showman et al. , 2009 ) 



822 We note that the condensation temperature of forsterite and enstatite is 



823 higher than 1,300 K at 1 mbar (Visscher et al. , 2010). Because silicon clouds 
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do not form, the NaoS clouds cannot form either. Further, the formation of 



[jNfa2S rehes on H2S, which is dissociated above the 1 mbar level (Zahnle et 



al.l 2009). This implies that any depletion of Na at high altitudes is most 



likely due to photoionization and/or thermal ionization. 

Figure [9] shows a dayside temperature (T-P) profile for HD209458b based 

on 



Showman et al. (2009). This profile is similar to the dayside temperature 



830 profile adopted by Moses et al. (2011 ). The figure also shows the condensation 



curves for forsterite and enstatite. The T-P profile crosses the condensation 
curve for forsterite below the 100 bar level. However, the temperature profile 
in the deep atmosphere is uncertain, and the current profile only barely 
crosses the condensation curve. Also, the formation of forsterite ties only a 



835 fraction of the total abundance of silicon into clouds (Visscher et al. , 2010). 



However, the T-P profile crosses the condensation curves for both forsterite 
and enstatite in a 'cold trap' near 10 mbar. To prevent this, the temperature 
in the cold trap would have to be T > 1,600 K. This is not totally unbelievable 
but probably unlikely. In any case, the temperature is close enough to the 
condensation curves so that moderate vertical transport might be able to 
preserve silicon above the cold trap. 



Spiegel et al. (2009) explored a range of turbulent diffusion coefficients 



Kzz that would be required to prevent the condensation of TiO and VO in the 
atmospheres of different EGPs, including HD209458b. In fact, they assumed 
that condensates form in the cold trap but are then transported to higher 
altitudes where they evaporate. Their work ignores the detailed chemistry 
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Figure 9: The dayside temperature-pressure (T-P) profile of 



HD209458b from Showman et al. (2009) that includes a temper- 



ature inversion at low pressures. The condensation curves for 
forsterite (dotted line) and enstatite (dashed line) are shown. 

847 of condensation, and thus the results are simply based on the ratio of Kzz to 

848 the diffusion coefficient estimated from 



849 where Vp is the particle settling velocity and H is the pressure scale height. 

850 The formation of condensates is probably too complicated for such a sim- 

851 plistic treatment, but the results provide some guidance on the value of Kzz 

852 that is required to lift the condensates from the cold trap. 

853 We calculated Dp for forsterite grains with a radius of 0.1 /im and densitjj^ 



^This is the density of the material in the particles, not the density of the particles in 
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854 of 3,200 kg m-3 ( |Fortney et alj |2003t [Cooper et alj |2003D . The settling 

855 velocity for such particles in the cold trap is fp ~ 3 x 10^^ m s^^ and 

856 Dp ^ 2 X 10^ s^^. Given that the upper edge of the cold trap is near the 

857 1 mbar level where Dp is higher, K^z ^ 10^ s~^ is sufficient to prevent the 

858 settling of the cloud particles in the lower atmosphere. We note that mass 

859 loss does not help to enhance the mixing of the particles at low altitudes 

860 significantly. The vertical velocity based on the mass loss rate of 10'^ kg s~^ 

861 is only 4.8 x 10^'' m s~^ at 10 mbar and 8 x 10^'^ m s^^ at 1 /ibar. 

862 Estimating K^z on extrasolar planets is very difficult, partly because there 

863 is no agreement on exactly what physical processes this parameter describes 

864 even in much more sophisticated solar system applications. The most recent 

865 estimates for HD209458b are based on assuming that Kzz ~ vH where v 



866 is the rms vertical velocity from circulation models (e.g.. Showman et al. 



2009). Based on the GCMs of Showman et al. (2009) and an assumed den 



868 sity dependence with altitude, Moses et al. (2011) estimated that the high 



869 pressure value of Kzz = 10^ m^ s~^, which implies that Kzz ^ 10^ m^ s~^ 

870 at p < 10 mbar. If such high values are realistic, turbulent mixing is almost 

871 certain to prevent the settling of the condensates and to preserve gaseous 

872 silicon in the upper atmosphere. 

873 We note that the above requirements on the value of Kzz may in fact be 

874 overestimated because they are based on the assumption that clouds particles 



the atmosphere. 
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875 form in the cold trap. Cloud formation has to be studied in the context of 

876 a photochemical model that includes the chemistry of condensation. If the 

877 chemical timescale is longer than the transport timescale, the cloud droplets 

878 may not form in the first place before the gas escapes from the cold trap. 

879 Also, the temperature structure near the cold trap needs to be constrained in 

880 greater detail. The formation of condensates is a complex problem that will 

881 be studied in future work (Lavvas et al., in preparation) in order to better 

882 constrain the required values of K^z- For our purposes it is sufficient to note 

883 that the current estimates of K^^ ™ply) in agreement with the observations, 

884 that sihcon clouds do not form in the upper atmosphere of HD209458b. 

885 4. Discussion and Conclusions 

886 We have used multiple observational constraints and theoretical models 

887 to characterize the upper atmosphere of HD209458b. Contrary to many of 

888 the earlier studies, we did not treat the thermosphere independently of the 

889 rest of the atmosphere. In fact, we have shown that observations of the upper 

890 atmosphere can be used to obtain useful constraints on the characteristics 

891 of the lower atmosphere. This is important because the extended upper 

892 atmospheres of close-in EGPs produce much larger transit depths than those 

893 arising from the lower atmosphere. In this work, we concentrated mostly on 

894 the FUV transit measurements. Future work should explore the possibility 

895 of extending the range of possible observations to other wavelength regions, 

896 as well as obtaining repeated observations in the FUV lines. Theoretical 

50 



models should be developed to support new observations and to clarify the 
interpretation of the existing data. 

In agreement with KIO, we showed that the H Lyman a transit observa- 
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2007, 


2008 



Hosseini, 2010) can be fitted with a layer of H in the thermosphere that is de- 
scribed by three simple parameters. These are the pressure at the bottom of 
the H layer, the mean temperature in the thermosphere, and a cutoff level due 
to ionization. The most important parameters are the pressure at the lower 
boundary and the mean temperature. Because H is the dominant species 
in the thermosphere, the data can be used to estimate the temperature of 
the thermosphere. Choosing a lower boundary pressure of 1 /ibar based on 
the location of the H2/H dissociation front in recent photochemical models 



(Moses et al. , 2011) and observational constraints (France et al. , 2010), we 
measured a mean temperature of about 8,250 K in the thermosphere below 3 
R.p. However, the uncertainty of the observations is large, and the la upper 
and lower limits on this temperature are approximately 6,000 K and 11,000 
K, respectiveljj^ 

We used a hydrodynamic model that treats the heating of the upper at- 
mosphere self-consistently and the average solar XUV spectrum (Paper I) to 
show that a mean temperature of 8,250 K in the upper atmosphere below 
3 i?p is higher than the maximum temperature allowed by stellar heating. 



"^This uncertainty does not include the possible uncertainties in the other parameters 
of the fit. 
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918 Given that a net heating efficiency equal to or higher than 100 % is unreahs- 

919 tic, this temperature requires either a non-radiative heat source, additional 

920 opacity, or it implies that the XUV flux of HD209458 is higher than the 

921 corresponding solar flux. Interestingly, this would also imply that the mass 

922 loss rate could be higher by a factor of 2 or more than previously anticipated 

923 (e.g.. Paper I). However, the uncertainty in the H Lyman a observations also 

924 allows for a slightly lower temperature of 7,200 K that is typical of stellar 

925 heating based on the average solar XUV flux and our best estimate of the 

926 net heating efficiency. Therefore the temperature implied by the H Lyman a 

927 observations and the mean temperature of the basic hydrodynamic models 

928 are in good agreement. 

929 We note that the temperature in the lower thermosphere near the 1 /ibar 

930 region has been constrained previously by Vidal-Madjar et al. ( 2011a|[b ) who 

931 used the Na D lines to constrain the density and temperature profiles in the 

932 atmosphere of HD209458b. Their results point to a temperature of ~3,600 

933 K that is actually higher than the temperature at the lower boundary of 

934 our hydrodynamic model and consistent with a high mean temperature at 

935 lower pressures in the thermosphere. However, we caution the reader that 

936 the temperature proffie based on the Na D lines may not be accurate. This is 



937 because Vidal-Madjar et al. (2011a) used a simple expression for the optical 

938 depth due to Na that is based on the scale height of the atmosphere [their 

939 equation (1)]. This expression is only valid if Na is uniformly mixed with 

940 H2. Since the authors also argue that Na is depleted (i.e., its mixing ratio 
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941 changes with altitude) above ~10 mbar, its density scale height cannot be 

942 used to estimate temperatures accurately. 

943 The detection of O in the thermosphere allowed us to constrain the mass 

944 loss rate from HD209458b based on the crossover mass concept formulated 



945 by Hunten et al. (1987). This is because O is transported to high altitudes 

946 in the upper atmosphere primarily by collisions with H. We found that a 

947 minimum mass loss rate of 6 x 10^ kg s~^ is required to prevent the diffusive 

948 separation of O. This is one of the most reliable constraints on the mass 

949 loss rate available at present. Our hydrodynamic calculation based on the 

950 average solar XUV flux predicts a mass loss rate close to 5 x 10'' kg s~^. 

951 This implies that species with a mass up to ~130 amu are uniformly mixed 

952 in the thermosphere. Similar constraints do not apply to heavy ions. They 

953 are transported to high altitudes by Coulomb collisions with that are 

954 much more efficient in preventing diffusive separation compared to collisions 

955 of neutral atoms with H. 

956 In agreement with KIO, we found that the presence of O with a solar abun- 

957 dance and a temperature based on the H Lyman a measurements explains 

958 the transits observed in the O I lines. Our models predict a line-integrated 

959 transit depth of approximately 4 % that deviates from the measured transit 

960 depth by 1.5cr and thus agrees with the uncertainty of the observations. How- 

961 ever, the predicted transit depths fall systematically short of the measured 

962 value. KIO suggested that a higher transit depth is possible if the 0/H ratio 

963 is supersolar, the temperature of the thermosphere is higher than expected, 
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964 and/or the observations probe escaping gas outside the Roche lobe of the 

965 planet, and we agree with their conclusions. As we explain below, we found 

966 that the same is true for the other heavy species. 

967 In order to calculate predicted transits in the C II and Si III lines, we 

968 created model emission line profiles for HD209458 based on SUMER ob- 



969 servations of the Sun (Curdt et al. , 2001) adjusted to the observations of 



970 HD209458 (Linsky et al. , 2010). We did not find evidence for significant ab- 

971 sorption by the ISM in the C II 1335.7 A or the Si III line. Parts of the C II 

972 1334.5 A line, on the other hand, are optically thick in the ISM and we took 

973 this into account in our models. We note that resolved observations of the 

974 emission line profiles can be used to characterize the composition of the ISM 

975 and the activity of the host star. This information is a valuable byproduct 

976 of the FUV transit observations that typically need to be repeated several 

977 times. 

978 With solar abundances, the same models that agree with the H Lyman a 

979 transit observations tend to underestimate the transit depths in the C II 

980 and Si III lines. Similarly with the O I lines, higher transit depths in the C 

981 II and Si III lines are possible if the mean temperature of the absorbers is 

982 higher than expected and/or the C/H and Si/H ratios are supersolar. With 

983 solar abundances, a \a agreement with the observed line-integrated transit 

984 depths is possible if the stellar XUV flux (or the stellar flux combined with 

985 an additional heat source) is higher than or equal to 2 times the average 

986 solar XUV flux. This corresponds to typical solar maximum conditions. 
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987 and it is quite interesting that a similar enhancement may be required to 

988 explain the relatively high mean temperature in the thermosphere that we 

989 estimated earlier. Alternatively, with heating based on the average solar 

990 flux, an agreement with the observations is possible with the 0/H, C/H and 



991 Si/H enhanced by a factor of ~3-5 relative to solar abundances (Lodders 



2003). In any case, the atmosphere is escaping with a minimum mass loss 



993 rate given above. This is evident because of both the high temperature of 

994 the thermosphere and the detection of heavy species at high altitudes. 

995 We note that the transit observations are affected by stellar variability 

996 that can lead to significant changes in the observed transit depths. Spatial 

997 variations of intensity on the stellar disk during maximum activity or limb 

998 brightening may render the transit occasionally undetectable, or enhance 

999 it by a significant factor. Generally, observations obtained during periods 



1000 of minimum activity are more reliable. Ben-Jaffel (2007) characterized the 

1001 short-term variability of HD209458 in the H Lyman a line, but the variabil- 

1002 ity of the O I, C II, and Si III are poorly characterized. This introduces 

1003 an additional element of uncertainty into the transit depths that needs to 

1004 be constrained by repeated observations of the star and the transits in the 

1005 FUV lines. It also means that the qualitative explanation of the present ob- 

1006 servations that is based on heating by the average solar XUV flux and solar 

1007 abundances (Paper I) may be sufficient even if the predicted transit depths 

1008 in the O I, C II, and Si III lines do not exactly match the current measure- 

1009 ments. On the other hand, if higher transit depths are confirmed, they can 



55 



1010 be used to further constrain the mean temperature and abundances as we 

1011 have shown. 

1012 The detection of heavy ions escaping the atmosphere can potentially be 

1013 used to constrain the magnetic field strength of the planet. We estimated an 

1014 upper limit of 0.04 mj for the magnetic moment that allows the heavy ions 

1015 to escape unimpeded at the equator. The estimated magnetic moment agrees 

1016 with the scaling laws for the magnetic field strengths of tidally locked close-in 



1017 EGPs by Griesmeier et al. (2004). On the other hand, a strong magnetic field 

1018 inhibits the flow of ions from the equator and may only allow for escape at the 



1019 poles (e.g., Trammell et al. , 2011). If the neutral-ion collision frequencies are 

1020 sufficient, the trapped ions may also prevent the neutral atoms from escaping. 

1021 We note that a uniform upward flux is required to mix the atmosphere, and 

1022 escape at the poles may not be sufficient to create a large enough obstacle to 

1023 explain the transits in the O I, C II, and Si III lines. Detailed models of the 

1024 magnetosphere that include the heavy species are required to assess if this is 

1025 the case or not. 

1026 The detection of Si^"*" in the upper atmosphere means that silicon can- 

1027 not condense to form enstatite, forsterite, or other condensates in the lower 

1028 atmosphere. This is clear because at least a solar abundance of silicon in 

1029 the thermosphere is required to explain the large transit depth in the Si III 

1030 line. According to the calculated temperature profiles for IID209458b (e.g.. 



Showman et al. , 2009), condensation is expected in a cold trap near the 10 



mbar level. Provided that the temperature is not much higher than expected 
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near the cold trap, efficient mixing is required to prevent condensation. Fol- 



1034 lowing an argument similar to that of Spiegel et al. (2009), we estimated 



that an eddy mixing coefficient of K^z ^ 10^ m^ s~^ below 1 mbar is suffi- 
cient to prevent the condensation of forsterite and enstatite in the cold trap. 
We note that much higher values than 10^ m^ s~^ were assumed by recent 



1038 photochemical models by Garcia Munoz (2007) and Moses et al. (2011) 



A stratospheric temperature inversion may also be necessary to suppress 
condensation. This is because the cold trap cannot extend to much lower 
pressure than 1 mbar or the required values of K^z become too high. Also, 
the lack of condensation implies that the temperature of the lower atmosphere 
should be relatively high. The detection of silicon in the upper atmosphere 
thus provides further evidence for a stratosphere on HD209458b that was 



1045 ffist proposed by Knutson et al. (2008). We note that existing radiative 



transfer models do not account for molecular opacity at UV wavelengths or 



visible absorbers potentially generated by photochemistry (e.g., Zahnle et al 



2009). Our results provide motivation for more detailed models of thermal 



structure below the thermosphere that can constrain the chemistry of the 
lower atmosphere. Once this is achieved, better constraints on the dynamics 
can be derived. 

We also address an old problem related to the atmosphere of HD209458b. 



1053 Based on the observed transits in the Na D lines, Charbonneau et al. (2002) 



argued that Na is depleted in the atmosphere because their cloudless so- 
lar composition model predicted significantly deeper absorption in the D 
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1056 lines. In addition to photoionization, molecular chemistry, and low primor- 

1057 dial abundance of Na, they suggested that the formation of high altitude 

1058 clouds can explain the observed depletion. Later Sing et al. ( 2008a|[b ) found 



1059 further evidence for the depletion of Na above the 3 mbar level, and argued 

1060 that condensation of Na2S is the most likely explanation. The fact that 

1061 silicon does not condense implies that condensation of Na2S is also unlikely. 

1062 The observed depletion is therefore most likely due to photoionization and/or 

1063 thermal ionization. However, the density profile and ionization state of Na 

1064 should be studied in the context of molecular and ion chemistry below the 

1065 0.1 /ibar level to verify that this is the case. 
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